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Abstract 

We propose a model equation for the dynamics of tree density in savannas. It only considers long-range competition 
among trees and the effect of fire, and despite short-range facilitation is not taken into account, we recover the 
standard spatial patterns observed in realistic situations. Therefore we show that facilitation is not a key ingredient 
in the modelling of savannas. We also study how fires affect the survival of trees in the long time limit as well as 
the maintenance of the patterns. 
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1. Introduction 

Savannas ecosystems are characterized by the long- 
term coexistence between a continuous grass layer and 
scattered or clustered trees [1]. Occurring in many 
regions of the world, in areas with very different cli- 
matic and ecological conditions, the spatial structure, 
persistence, and resilience of savannas have long in- 
trigued ecologists . However, despite substantial 
research, the origin and nature of savannas have not 
yet been well explained. 

Savanna tree populations often exhibit pronounced, 
non-random spatial structure 1151- 41 Cil . Much research 
has therefore focused on explaining how spatial pat- 



terning in savannas arises ||7|-|yi [1JJ, [12|]. It has fre- 
quently been claimed that spatial pattern in savanna 
vegetation can be explained by a combination of long- 
distance inhibition and short-distance facilitation II 1 3l — 
This combination of mechanisms is known to 
produce patterns in many other natural systems. While 
the facilitation component has often been thought to 
be a key component in previous savanna studies |9l 
3 [13, El , Rietkerk and Van de Koppel H, based on 
a review of spatial patterns observed in diverse ecosys- 
tems, suggested that only the long-distance competi- 
tion part is key to pattern formation. To our knowl- 
edge, this explanation has not been explored for a sa- 
vanna system. On the other hand, a similar mechanism 
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of competitive interactions between species, similar to 
the one presented here for savannas, has been shown 
to give rise to clusters of species in the context of clas- 
sical ecological niche theory. Scheffer and van Nes 
lfl9ll showed that species distribution in niche space 
was clustered, and Pigolotti et al. 12011 showed that 
this arises as an instability of the nonlocal nonlinear 
equation describing the competition of species. 

In the study of savannas, exploring whether or not 
long-range competition, acting nonlinearly alone or in 
combination with some others mechanisms, can ac- 
count for spatial patterning is important because the 
issue of short-range facilitation among savanna trees 
has been highly contentious. Some authors have pre- 
sented evidence of short-distance facilitation @], 
while others have demonstrated evidence of short- 
range competition Jj, |7l [10J], with conflicting reports 
sometimes arriving from the same regions. Indeed, 
the most frequently invoked mechanism of short-range 
facilitation in savannas, that tree clumps promote 
infiltration of rainwater that would instead become 
runoff on patches of bare or sparsely- vegetated ground 
1 15t|21[], seems unlikely to apply to all savannas that 
display spatial pattern. The local infiltration mech- 
anism seems most plausible for extremely arid sys- 
tems, where there is substantial bare ground that acts 
as a barrier to infiltration. For mesic to wet savan- 
nas, the infiltration mechanism may not be as impor- 
tant because much more of the ground is covered by 
vegetation. Under these conditions, however, substan- 
tial competition among trees may still occur for water 
or other limiting resources. The root architecture of 
savanna trees and shrubs, with lateral roots of many 
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species extending several meters in each direction, 
provides a clear physical mechanism by which non- 
local competition for water or other resources could 
be mediated iflH. 

Here we develop a minimalistic model of savannas 
that considers two of the factors that are thought to 
be crucial to structure semi-arid to mesic savannas: 
tree-tree competition and fire. Our main focus is on 
spatially nonlocal competition, and employing stan- 
dard tools used in the study of pattern formation phe- 
nomena in physics (stability analysis and the structure 
function), we explore the conditions under which it 
can produce non-homogeneous spatial distributions. 
We also include the effect of fires as a function of 
the local density of trees. We will show how impor- 
tant fires are in the survival of vegetation, i.e., in the 
long-term coexistence of trees and grass. On the other 
hand, fire does not play a crucial role in the emer- 
gence of spatial structures in our model: we observe 
the same sequence of patterns when competition in- 
tensity is manipulated, independently of the frequency 
of fires. Competition is therefore the main mechanism 
governing the formation of spatial structures. We fur- 
ther examine the role of demographic stochasticity in 
modifying both spatial patterns and the conditions un- 
der which trees persist in the system in the presence 
of fire. Our model is inspired by the one presented by 
Calabrese et al. in II ill , and complements their work, 
providing further analytical results. 

Our results clearly demonstrate that this simple sys- 
tem, with no short-range facilitation, can produce the 
full spectrum of spatial patterns reported from mod- 
els employing both short-range facilitation and long- 
range inhibition (competition). We conclude by dis- 
cussing the implications of these results for the debate 
on whether the balance of short-range processes af- 
fecting savanna trees is positive, negative, or is vari- 
able among systems. To complete our work we in- 
clude an appendix where we study the effect of exter- 
nal fluctuations (mimicking for e.g. rainfall) on sa- 
vanna dynamics. 



2. The Deterministic Model 



In this section we derive the deterministic evolution 
equation for the local density of trees, such that dy- 
namics is of logistic type and we only consider tree- 
tree competition and fire. We study the formation of 
patterns via stability analysis and provide numerical 
simulations of our model, showing the emergence of 
spatial structures. 



2.1. The nonlocal savanna model 



In 111 111 it is introduced a simple discrete particle 
model on a lattice for a savanna system that considers 
the birth-death dynamics of trees, and where tree-tree 
competition and fire are the principal ingredients. It is 
derived a deterministic differential equation describ- 
ing the time evolution of the global density of trees, 
p(f), that reads: 



^ =bP E (p)p(t)(l-p(t))-ap(t), 
at 



(1) 



where a and b are the death and birth rates, respec- 
tively, and Pe(p) is the probability of establishment 
given that a seed has been dispersed. Because compe- 
tition and fire are independent, P E - PcPf, where Pq 
is the probability of surviving the competition, and Pf 
is the probability of surviving a fire event. 

Generalizing Eq. ([TJ, we propose an evolution equa- 
tion for the space-dependent (local) density of trees, 
p(x, t): 



dp(x, t) 
dt 



= bP E p(x, 0(1 - p(x, f)) - ap(x, t). (2) 



We allow the probability of overcoming competition 
to depend on tree crowding in a local neighborhood, 
decaying exponentially with the density of surround- 
ing trees as 



Pc = exp 



G(x - r)p(r, t)dr\ 



(3) 



where 6 is a parameter that modulates the strength of 
the competition, and G(x) is a positive kernel function 
that introduces a finite range of influence. Pf depends 
on the local density of grass, 1 - p(x, t), as 



(4) 



cr + 1 - p(x, t) 



where cr governs the resistance to fire. Notice how our 
model is close to the one in [11] through the defini- 
tions of Pc and Pf. The deterministic evolution equa- 
tion that considers tree-tree competition and fire for 
the spatial tree density is 



<5p(x, f) 
dt 

where 



b eff (p)p(x, t) (1 - p(x, f)) - ap(x, t), (5) 



beffip) = 



(6) 



<T + 1 - p(x, t) 

Thus, we have a logistic-type equation with an ef- 
fective growth rate that depends nonlocally on the den- 
sity itself. 

One could include a diffusion term like D(p)V 2 p in 
Eq. (|5]l to account for local dispersal of trees (that may 



2 



reflect local seed dispersal giving rise to adult trees) 
123 , 24]. Though we are interested in isolating the ef- 
fect of competition among trees and fires, its possible 
consequences are discussed in the section[4] 

As noticed before, the idea behind the nonlocal 
competition is to capture the effect of the long roots 
of a tree. Thus the nonlocal interaction scale is of the 
order of the length of the roots [2]. Since the roots are 
the responsible for the adsorption of resources (wa- 
ter and soil nutrients), a strong long-range competition 
term implies strong resource depletion. 

Since we are interested in the effect of competition 
and fire on the distribution of savanna vegetation, we 
will try to fix all the parameters but cr and 6. Firstly, 
(as in [11]) the death rate is taken a = 1, which can be 
accomplished by properly rescaling time. We can also 
fix the growth rate b by considering the homogenous 
solution of Eq. (0, po = po{b, cr, 6). It can be obtained 
by solving 



b eff (p )(l -po)- 1 =0, 



(7) 



with real data from African savannas featuring in- 
termediate to high mean annual precipitation levels 
1I3 H25I1 . for which competition for water is no longer 
limiting (6 = 0) and fire is not an issue (cr — > 00). 
In this limit, the upper tree cover has been measured 
around po = 0.8, and the stationary solution of the 
model depends only on the birth rate 



Po(o- ■ 



3) 



b- 1 



(8) 



which yields b — 5 solving for a fixed po = 0.8 ( Hill ). 
In the following we just consider the dependence of 
our results on S and cr. In particular, po = po(cr, 6). 

One of the characteristics of this model that makes it 
suitable for modelling savannas is the long-term coex- 
istence between grass and trees, characteristic of these 
ecosystems and shown by our equation. To explore it, 
we study the long-time behavior of the system and an- 
alyze the stationary solutions of Eq. (0, which has two 
fixed points. The first one is the absorbing state rep- 
resenting the absence of trees, po = 0, and the other 
can be obtained, in the general case, solving Eq. (O 
numerically. In the regime where po is small (near the 
critical point), if competition intensity, 5, is also small, 
it is possible to obtain an analytical expression for the 
critical value of the probability of surviving a fire, cr*, 
which is independent of the strength of competition 



1 



b-1 



(9) 



i.e. cr* = 0.25 for b — 5. Outside of the limit where 
5 <K 1, we can solve Eq. (|7]l numerically in po to show 
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Figure 1 : Grass-coexistence phase transition. The lines come from 
the mean field solution, Eq. {7), and the dots from the numerical 
integration of the PDE, Eq. {3) in a region of 1 ha. Notice the finite- 
size effect in the shift of the critical point in the dots. Nevertheless, 
its value remains independent of S, as expected. 



that the critical value of the fire resistance parameter, 
cr*, also does not depend on competition. 

For higher fire survival probability a steady state 
with trees is stable (Fig. [TJ. The model, then, shows a 
transition from a state where grass is the only form of 
vegetation to another state where trees and grass coex- 
ist at cr*. 

2.2. Linear stability analysis 

The spatial patterns appearing in the nonlocal sa- 
vanna model can be studied by performing a linear 
stability analysis J26ll of the stationary homogeneous 
solutions of equation ©, po = po(cr, 6). The stability 
analysis of po is performed by considering small har- 
monic perturbations around it, p(x, f) — Po + ee^ - *'", 
e <k 1 . After some calculations ( Appendix A 1, one 
arrives at the dispersion relation 



A(k;cr,6) = b eff (p ) 1 -^^ 



-beff(po) 



po[2-po+SG(k)(p - 1 )jpp - 1 -a)] 
(cr-po+1) 



- 1, 



(10) 



where G(k), k 
kernel, 



Ikl, is the Fourier transform of the 



G<10 



G(x)e ikx dx. 



(11) 



The critical values of the parameters of the tran- 
sition to pattern, 5 C and cr e , and the fastest growing 
wavenumber k c are obtained from the simultaneous 
solution of 



A(k c ; cr c , <y = 0, 
= 0. 



(12) 
(13) 



Note that k c represents the most unstable mode of 
the system, which means it grows faster than the oth- 
ers and eventually dominates the state of the system. 
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Therefore, it determines the length scale of the spa- 
tial pattern. These two equations yield the values of 
the parameters 6 and cr at which the maximum of the 
curve A(k), right at k c , starts becoming positive. This 
signals the formation of patterns in the solutions of 
Eq. ©. As Eq. (TT~3T > is explicitly written as 

A'(k c ) = b eff (p )6poG'(k c )(p - 1), (14) 

the most unstable wavenumber k c can be obtained by 
evaluating the zeros of the Fourier transform of the 
kernel. 

Equation (TlOb shows that competition, through the 
kernel function, fully determines the formation of pat- 
terns in the system. If the Fourier transform of G never 
takes positive values then A(k; cr, 6) is always negative 
and only the homogenous solution is stable. However, 
when G can take negative solutions then patterns may 
appear in the system. What does this mean in biolog- 
ical terms? Imagine that we have a family of kernels 
described by a parameter p: G(x) = exp(-|(x) / R\ p ) {R 
gives the range of competition). The kernels are more 
peaked around x = for p < 2 and more box-like 
when p > 2. It turns out that this family of functions 
has non-negative Fourier transform for < p < 2, 
so that no patterns appear in this case. A lengthy dis- 
cussion of this property in the context of competition 
of species can be found in M20I1 . Thus, the shape of 
the competition kernel dictates whether or not patterns 
will appear in the system. If pattern formation is pos- 
sible, then the values of the fire and competition pa- 
rameters govern the type of solution (see below). 

Our central result for nonlocal competition is that, 
contrary to conventional wisdom, it can, by itself pro- 
mote the clustering of trees. Whether or not this occurs 
depends entirely on the shape of the competition ker- 
nel. For large p we have a box-like shape, and in these 
cases trees compete strongly with other trees, roughly 
within a distance R from their position. The mech- 
anism behind this counterintuitive result is that trees 
farther than R away from a resident tree area are not 
able to invade it (their seeds do not establish there), so 
that an exclusion zone develops around it. For smaller 
p there is less competition and the exclusion zones dis- 
appear. 

For a more detailed analysis, one must choose an 
explicit form for the kernel function. Our choice is 
determined by the original Pq taken in lUHl . so that 
it decays exponentially with the number of trees in a 
neighborhood of radius R around a given tree. Thus, 
for G we take the step function (limit p — > oo) 



where R is the interaction range that determines the 
competition distance among trees. For this kernel the 
Fourier transform is R27I1 G(k) = 27iR 2 J\{kR)lkR and 
its derivative is G'(k) = -nR 2 J2(kR)/k, where k = |k|, 
and Jj is the /'"-order Bessel function. Since G(k) 
can take positive and negative values, pattern solutions 
may arise in the system, that will in turn depend on the 
values of 5 and cr. The most unstable mode is numer- 
ically obtained as k c = 5.136//?. Because a pattern of 
n cells is characterized by a wavenumber k c = 2nnjL, 
where L is the system size, the typical distance be- 
tween clusters d, - Lin is given by d, w 1.22/?, i.e., it 
is approximately the range of interaction R. This result 
is independent of the other parameters of the system. 

The phase diagram of the model, computed numer- 
ically, is shown in Fig. [2] where depending on 5 and cr 
we plot the spatial character of the steady solution (ho- 
mogeneous or inhomogeneous). Note that increasing 
competition enhances the inhomogeneous or pattern 
solution. This is because, as we are now in the case of 
a kernel giving rise to clusters, increasing 6 makes it 
more difficult to enter the exclusion zones in-between 
the clusters. For very strong competition (high, unre- 
alistic, 5), fire has no influence on the pattern. How- 
ever, for small to medium values of 6, frequent fires 
(i.e., small cr) enhances clustering since it also makes 
more difficult to enter into exclusion areas. 

The critical line separating these two solutions (pat- 
tern and homogeneous) can be analytically computed 
as a function of the parameters 6, cr, po and G(k c ) (see 
Eq. ( IB. II ) in Appendix Bi. In Figure [2] we have plot- 
ted (with crosses) this critical line separating homoge- 
neous and pattern solutions for the step kernel. Note 
that the stationary density of trees, po, must be com- 
puted numerically from Eq. (0. 

With b — 5, in the absence of fire (cr — > oo), and for 
weak competition, we can take the limits 5 — > and 
cr — > oo of the dispersion relation Eq. ( [Tol l, leading to 



A(k; ci->0,cr-^oo)=4 - 10p , 



(16) 



G(|r|) 



1 if |r| < R 
if \t\>R, 



(15) 



which is negative for the typical density of trees in the 
system for these values of the parameters (See Fig.[TJ). 
This means that in this limit, trees are uniformly dis- 
tributed in the system as there is no competition, and 
space does not play a relevant role in the establish- 
ment of new trees. This result agrees with the phase 
plane plotted in Figure [2] In biological terms, there 
are no exclusion zones in the system because there is 
no competition. 

2.3. Numerical simulations 

The previous analysis provides information, de- 
pending on the competition and fire parameters, about 



4 



10 



HOMOGENEOUS 

p o = unstable 

* 




HOMOGENEOUS 5 

p = stable 

Figure 2: Possible solutions of the mean field equation Eq. {3) for 
b = 5.0, a = 1.0, and a step kernel. The absorbing-active transition 
in the mean field model is shown at <x* with circles (o), as well as 
the homogeneous-pattern transition, indicated with crosses (x). The 
diamond, the square, and the up-triangle show the value of the pa- 
rameters cr and S taken in Figures l3]4l and[5]respecfively. The stars 
point out the transition to inhomogeneous solutions in the stochastic 
model as described in Section|3] with T = 0.2. 



when the solution is spatially homogenous and when 
trees arrange in clusters. However, the different shapes 
of the patterns have to be studied via numerical sim- 
ulations lr28ll of the whole equation of the model. We 
have taken a finite square region of savanna with an 
area of 1 ha., allowed competition to occur in a circular 
area of radius R — 8 m, and employed periodic bound- 
ary conditions and a finite differences algorithm to ob- 
tain the numerical solution. Similarly to what has been 
observed in previous studies [15, 11711 . different struc- 
tures: gaps, stripes, and tree spots, are obtained in the 
stationary state as we make competition stronger given 
a fixed value of the fire frequency or, on the other hand, 
as we decrease the resistance to fires for a given com- 
petition intensity. In both equivalent cases, we observe 
this spectrum of patterns as far as we go to a more dry 
state of the system, where resources (mainly water) are 
more limited (see Figs. OH). It is notable that our mod- 
eling considers only nonlocal competition, and not, as 
already commented, short-range facilitation and long- 
range competition mechanisms as is typical in other 
studies ni5[|l7ll . which gives an idea of the robustness 
and relevance of these type of spatial patterns. 

This same sequence of appearance of patterns has 
been already observed in [ 15] and indicates that, for a 
given value of the fire parameter, increasing 5, i.e. de- 
creasing the probability of surviving competition, new 
trees cannot establish into the exclusion areas so clus- 
tering is enhanced. Self-organized vegetation patterns 
can persist far into regions of high aridity where plants 
would become extinct if homogeneously distributed 
115]. 

A much more quantitative analysis of the periodic- 




Figure 6: Structure function defined in equation )17t for different 
values of the demographic noise intensity. The main peak in A = 10 
reveals the periodicity of the solution and the equivalence among 
patterns independent of the fluctuations. 5 = 9.8, a = 2.9, R = 8, 
a= 1.0,6 = 5.0. 



ity in the patterns can be performed via the structure 
function. This will be helpful to check the previous re- 
sults and, especially, for the analysis of the data of the 
stochastic model of the next section, for which we will 
not present analytical results. The structure function is 
defined as the modulus of the spatial Fourier transform 
of the density of trees in the stationary state, 



S(k) 



dxe p(x, t -» 00) 



(17) 



where the average is a spherical average over the 
wavevectors with modulus k. The structure function 
is helpful to study spatial periodicities in the system, 
similar to the power spectrum of a temporal signal. 
Its maximum identifies dominant periodicities, which 
in our case are the distances between tree clusters. 
Note that the geometry of the different patterns can- 
not be uncovered with the structure function, since 
this involves a spherical average. For our determin- 
istic savanna model, the peak is always at A c = 10m 
(Figure [6] with F = 0; this parameter will be defined 
later on), which is in good agreement with the theo- 
retical result provided by the linear stability analysis 
A — 2njk max — 9.78m. 

3. Stochastic model 

The perfectly periodic patterns emerging from the 
current deterministic model seem to be far from 
the disordered ones usually observed in aerial pho- 
tographs and shown by individual based models ll ill . 
Thus, to provide a reasonable description of the un- 
derlying small-scale birth and death dynamics, we add 
demographic noise to the standard deterministic equa- 
tion, which will take into account the intrinsic fluctu- 
ations always present at the individual level. 
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Figure 3: Grass spot pattern in the determinis- 
tic model pointed out with a diamond in Fig- 
ure [2] The parameters are r5 = 7.0, cr = 2.9, 
R = 8.0, b = 5.0 and a = 1.0. 



Figure 4: Striped pattern in the deterministic 
model pointed out with a square in Figure [2] 
The parameters are 8 = 8.0, cr = 2.9, if = 8.0, 
2> = 5.0anda= 1.0. 



Figure 5: Tree spot pattern in the determinis- 
tic model pointed out with a triangle in Figure 
[2] The parameters are S = 11.0, cr = 2.9, 
R = 8.0, b = 5.0 and a = 1.0. 




Figure 7: Maximum of the structure function for different values 
of the competition parameter <5 at long times. The fire parameter 
is fixed at cr = 2.9. Black circles show the transition to pattern in 
the deterministic model while red squares refer to the model with 
demographic noise, T = 0.20. 



The interpretation of the field p(x, f) is the density 
of tree (active) sites in a small volume, V. If we 
think of trees as reacting particles which are born and 
die probabilistically, the number of reactions taking 
place in a volume, V, is proportional to the number 
of particles therein, TV, with small deviations. If TV is 
large enough, the central limit theorem applies to the 
sum of /V independent random variables and predicts 
that the amplitude of the deviation is of the order of 
V^V oc -y/p(x, f) 12911 . These fluctuations are referred 
to as demographic noise. The macroscopic equation is 
now stochastic, 



= b eff (p)[p(x,t)-p 2 (x,t)]- 



ap{x, t) + T Vp(x, t)t)(x, t), (18) 



dp(x, t) 
dt 



where T oc Jb e ff (but we take it as a constant, 1 30] ) 
modulates the intensity of 77(x, ?), a Gaussian white 
noise term with zero mean and correlations 
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Figure 8: Active-absorbing phase transition both in the deter- 
ministic model (Circles) and considering the demographic noise 
(Squares). In the case of the stochastic model, we integrate the 
Eq. QU with T = 0.2 and average the density of trees in the steady 
state. 



We first investigate the effect of intrinsic fluctua- 
tions on the persistence of trees in the system. We 
show (Fig. [8]) that the critical point, cr*, depends on 
the value of the competition parameter 5. Increasing 
the competition strength, which implies a more arid 
savanna, makes trees less resistant fires. As is seen 
in the figure this effect is rather small, so that when 6 
increases the transition to the grassland state appears 
only for a slightly larger cr (i.e., less frequent fire). 
The reason seems to be that fire frequency and inten- 
sity depend on grass biomass. Seasonally wet savan- 
nas support much more grass biomass that serves as 
a fuel for fires during dry season. Dry savannas have 
much lower grass biomass, so they don't burn as often 
or as intensely. 

We numerically explore the stochastic savanna 
model using an algorithm developed in J30ll . Note 
that the noise makes the transition smoother so the 
change from homogeneous to inhomogeneous spatial 
patterns is not as clear as it is in the limit where the de- 
mographic noise vanishes (See Fig. [7j. Figure [7] also 
shows how the maximum of the structure function de- 
cays with the competition strength when demographic 



fluctuations are considered. As we are keeping the re- 
sistance to fires fixed, increasing the competition in- 
tensity takes the system towards the absorbing phase, 
where the solution is homogeneous and the structure 
function shows no peaks. 

We have studied the dynamics of the system for 
some values of the fire and competition parameters, 
focusing on the role of the noise strength F. Despite 
the obvious increasing randomness in the patterns with 
T (Figures 19110b . they are statistically equivalent to 
those in the deterministic model (See the structure 
function in Figure |6j peaked around A = 10). There- 
fore, the model presents an active-absorbing transition 
with the noise strength, F, being the control parame- 
ter, and we obtain (more realistic) irregular patterns of 
trees in the whole clustered region, where the shape of 
the pattern in the deterministic model is not relevant. 
We also find (shown with stars in Fig. |2]l that the ap- 
pearance of patterns in the system now occurs at lower 
values of the competition parameter, 5, than in the no 
noise limit, but the shape of the transition line remains 
qualitatively similar. 



4. Discussion 

The striking observed spatial arrangements of sa- 
vanna trees have led to speculation on the mecha- 
nisms producing these patterns H OH 15]. A 
key step in such an analysis is defining the most par- 
simonious combination of mechanisms that will pro- 
duce the pattern in question. The literature on savanna 
tree spatial patterning with kernel-based models has 
tended to focus on the interplay between local facili- 
tation and long-range competition, under the assump- 



tion that both are essential for pattern formation ||F3 
[l^l . Following Rietkerk and van de Koppel's JH] ob- 
servation that not all of the ecosystems in their review 
of spatial patterning featured a clear short-range facil- 
itation mechanism, we have studied the conditions un- 
der which long-range competition alone could account 
for savanna tree spatial pattern. 

We have shown that nonlocal competition, acting 
without facilitation, is indeed sufficient to induce the 
full range of observed spatial patterns, provided the 
competition term enters nonlinearly in the equation 
for the density of trees, and that competition is strong 
enough, so resources in the system are limited. The 
key technical requirement for this effect to occur is that 
the Fourier transform of the competition kernel be- 
comes negative for certain values of the wavenumber, 
k. In practice this means that kernels whose Fourier 
transform takes negative values for some wavenumber 
values, will lead to competition driven clustering. The 



other mechanism we have considered for a minimalis- 
tic but realistic savanna model, fire, has been shown to 
be relevant for the coexistence of trees and grass and 
for the shape of the patterns. However, competition is 
the main ingredient allowing pattern solutions to exist 
in the model. If the shape of the kernel allows these 
types of solutions, then the specific values of fire and 
competition parameters determine the kind of spatial 
structure that develops. It is also worth to mention 
that one can observe the full spectrum of patterns in 
the limit where fires vanish (cr — > oo) as we increase 
the competition, but not when there is not competition, 
(5 = 0, regardless of the value of the fire term. There- 
fore, we conclude that the nonlocal competition term 
is the responsible of the emergence of clustered distri- 
butions of trees in the model, with the fire term playing 
a relevant role when it comes to fixing the value of the 
competition parameter at which patterns appear. 

These findings have the potential to shift the fo- 
cus of the debate over the role of short-range factors 
on large-scale savanna tree spatial pattern. Our re- 
sults demonstrate that short-range facilitation is not a 
necessary condition for pattern formation, suggesting 
that the precise balance of positive and negative short- 
range forces may have little or no effect on large-scale 
pattern in real savannas. While understanding the net 
effects of these short-range interactions among trees 
may be relevant for other aspects of savanna tree ecol- 
ogy, our results suggest that, when spatial pattern is 
the focus, emphasis should instead shift to identifying 
when, where, and to what extent nonlocal competition 
will be important. 

The other important feature of savannas, the charac- 
teristic long-time coexistence of trees and grass is well 
captured with our model (Figures [T] and [HJ. Besides, 
introducing demographic fluctuations, as it is shown 
in Figure [8] makes our approach much more realistic, 
since the persistence of trees in the face of fires is re- 
lated to the water in the system. On the other hand, 
demographic sctochasticity causes tree extinction at 
lower fire frequencies than in the deterministic case, 
because random fluctuations in tree density are of suf- 
ficient magnitude that tree density can hit zero even 
if the deterministic stationary tree density (for a given 
fire frequency) is greater than zero. 

5. Summary 

We have shown the formation of patterns in a min- 
imal spatially nonlocal savanna model, as well as the 
transition from a state where trees and grass coexist 
to an absorbing one, where only grass is present. The 
salient feature of the model is that it only considers 
nonlocal (and nonlinear) competition, which differs 
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Figure 9: Irregular pattern in the model with demographic noise. 
The parameters are r = 0.10, S = 11.0, cr = 2.9, R = 8.0, b = 5.0 
and a = 1.0. 



Figure 10: Irregular pattern in the model with demographic noise. 
The parameters are T = 0.20, S = 1 1.0, cr = 2.9, R = 8.0, b = 5.0 
and a = 1.0. 



from standard kernel-based savanna models that fea- 
ture both short-range facilitation and long-range com- 
petition. The same sequence of spatial patterns ap- 
pears in both approaches, confirming Rietkerk and van 
de Koppel's [8] suggestion that short-range facilitation 
may not be a necessary condition for spatial pattern 
formation. 



Inspired by 111 111 , we have proposed a nonlocal de- 
terministic macroscopic equation for the evolution of 
the local density of trees where fire and tree-tree com- 
petition are the dominant mechanisms. It is a nonlocal 
interaction model like the kernel-based models J2] for 
vegetation pattern formation in arid ecosystems, but 
while the kernel function is usually used to account 
for the short-range facilitation and long-range compe- 
tition, we use it exclusively to define the competition 
neighborhood (the length scale of the interaction). In 
other words, our approach only considers competitive 
interactions, without facilitation. If the kernel func- 
tion falls off with distance very quickly (the Fourier 
transform is always positive) the system only has ho- 
mogenous solutions. In the opposite case, patterns 
may appear depending on system parameters (fire and 
interaction strength), and in a sequence similar to the 
spatial structures appearing in standard kernel-based 
models. Under less favorable environmental condi- 
tions, trees tend to arrange in more robust structures 
to survive (Fig. |5). Biologically, trees are lumped in 
dense groups, separated by empty regions. Entrance of 
new trees in these exclusion zones is impossible due to 
the intense competition they experience there. 

This deterministic treatment is formal, and we have 
shown the different spatial distributions of the trees 
that occur as competition becomes more intense, con- 
cluding that self organization of trees is a good mech- 
anism to promote tree survival in extremely arid con- 
ditions 11511 . Trees tend to cluster in the high competi- 
tion (few resources) limit (Fig. [5]), due to the formation 
of exclusion zones caused by non-local competition, 



and not as a result of local facilitation. However, be- 
cause we are dealing with a deterministic model, the 
patterns are too regular and the transition between the 
grass-only and a tree-populated states is independent 
of tree competition. Fluctuations, due to the stochas- 
tic nature of individual birth and death events, are in- 
cluded for a more realistic description. Calabrese et 
al. Hill also noted that savanna-to-grassland transition 
was independent of competition intensity in the de- 
terministic case, but not in the stochastic case. In the 
present model, both the grassland to savanna transition 
and the spatial structures that develop are influenced 
by demographic stochasticity. 

We have quantified the characteristic spacing of 
spatial patterns through the structure function. The 
irregular patterns produced by the stochastic model 
still have a dominant wavelength whose value is the 
same as in the deterministic model and depends only 
on the value of the range of the interaction, R in the 
kernel function. The match between the typical spa- 
tial scale of the patterns and the characteristic distance 
over which nonlocal competition acts indicates that 
competition is responsible for the presence of clus- 
tered spatial structures. 
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Appendix A. Linear stability analysis 

This appendix shows the details of the linear stabil- 
ity analysis, in particular how the dispersion relation 
in Eq. (TT~0t is obtained. We consider the stationary so- 
lution po plus a small perturbation, 



p(x, f) = po + eif/(x, t), 



(A.l) 
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where e <K 1. Substituting Eq. dA.U into the orig- 
inal equation (0 and retaining only linear terms in 
e, and then taking harmonic perturbations of the type 
tf/(x, t) = exp(/if - ;k • x), we arrive to the relation dis- 
persion 



A(k) = bCcr{po - pi) 



i 



6{k)5 



(cr+l-po) 2 cr+l-po 



+bC<r-^- - a, 

cr+l-po ' 



(A.2) 



where G is the Fourier transform of the kernel and 
C = exp^-dpo J G(x-x')c/xJ. This expression can 
be written as Eq. ( TTOb using the definition of b e ff(po) 
and taking a - I. 



Appendix B. Expression of the transition to pat- 
tern critical line. 

We show here the analytical expression for the crit- 
ical line in the transition from homogeneous to in- 
homogeneous solutions. Starting from Eq. ( [T2l it is 
possible to write an expression for the value of the 
resistance to fires parameter, cr, at which the macro- 
scopic equation © starts showing pattern solutions, 
as a function of the competition parameter, 6, and the 
most unstable mode k c . Considering the value of the 
parameters taken in our study, b — 5 and a — 1, it is 

„ _ (po - l)[5(p - 1 )(SC(k c )p - 1 )-2e M2 >'Q ] 



" ° lo[l-2p„+<5G(i c )po(l+po)-e*' s2 ™/5] 
(Po-1) ^5[5(p -l)HSG(k c )p -iy--4e s '' R2 POp ] 
1 0[ 1 -2p +5G(k c )p a ( 1 +p )-e 4 " R2 "o /5] 



(B.l) 



This complicated expression must be evaluated nu- 
merically together with the solution of Eq. (|7]i for the 
stationary density of trees, which is also a function of 
the competition and fire parameters. We show the re- 
sults in Figure|2] where the curve, represented with the 
black crosses, fits perfectly with the numerical results 
from the linear stability analysis. 

Appendix C. The effect of rainfall: Random 
switching between death and birth 

One of the key ingredients for the long coexis- 
tence between grass and trees is the largely inhomoge- 
neous temporal distribution of precipitations over time 
|3l|28l]. We have studied this environmental variability 
following the idea in 01811 . considering the switching 
between unstressed vegetation growth, given by the 
first term in ©, and drought-induced vegetation de- 
cay, represented with the second term in Eq. ©. These 



processes take place each time step with probability P 
and 1 — P, respectively. From now on, we call 

f b [p(X, ?)] = b eff {p) [p(x, t ) - p 2 (x, ?)] , 

f d [p(x,t)] = -ap(x,t), (C.l) 

and 

/ ± bo(x, 0] = \ [fb[p(*, t)] ± f d [p(x, f)]] ■ (C.2) 

The random dynamics of the system is written in 
terms of a stochastic partial differential equation, 



dp(x, t) 

at 



= /+ [p(x, t)] + /_ [p(x, tm n (t), (C.3) 



where £/„(f) is a dichotomous noise (DMN), assum- 
ing values +1 (wet season) and -1 (dry season) with 
probability P and 1 — P, respectively. 

If the rate of random switching, here taken as the 
inverse of the integration time step, is relatively fast 
respect to the rate of convergence to equilibrium in 
each of the two states, we can replace the noise term 
in Eq. ( IC 3b with its average value, < £/„(f) >= 1 —2P. 
It is meaningful since the rainfall seasons are much 
shorter than the time needed to reach one of the equlib- 
rium stationary states of death and birth processes, 
p(x, t) — 0, 1, respectively. This substitution leads to a 
deterministic equation 



dp(x, t) 
dt 



f + \p(x,t)]+f-\p(x,t)](l-2P), (C.4) 



where we will be able to perform linear stability anal- 
ysis as usual. The new dispersion relation is easily 
obtained, 



A(k;cr,S)=b eff (po)P 1 -^^ 



(l-P) 



1 1 x n Po[2-po+i5G(*)(Po-l)(Po-l-°")+2o-] 

-b eff (po)P— (S^p7+T) > ( Q5 ) 

which means that the main effect of the dichotomous 
noise is to renormalize the rates a and b. The patterns 
observed now are the same as the ones in the deter- 
ministic case, though the regions where they emerge 
change in accordance with this renormalization. Thus, 
the effect of stochastic precipitation, as modeled with 
this random switching mechanism, is a change of the 
parameter values for the different transitions observed 
in the deterministic continuum model Eq. (0. 

According to the value of P, an absorbing-active 
phase transition is observed, P c m 0.20. Small values 
of P, meaning long dry season, lead to an absorbing 
state while increasing the probability of raining im- 
plies the appearence of trees in the system. In this 
latter case, the solution can be either homogeneous or 
showing spatial patterns, depending on fire and com- 
petition. 
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This attempt to model rainfall has not been very 
succesful and does not give a lot of new information. 
Much effort of future research should be put on this 
point, trying to get much more realistic modelling of 
external environmental variability, according to empir- 
ical observations, with long runs of dry years and rare 
wet years. 
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